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The influence of dispersion-less quantum optical phonons on the 
phase diagram of a quarter-filled Hubbard chain is studied using the 
Density-matrix renormalization group technique. The ground state 
phase diagram is obtained for frequencies corresponding to the intra- 
molecular vibrations in organic conductors. For high vibrational modes 
the system is only slightly affected by the electron-phonons coupling. It 
remains in a Liittinger liquid phase as long as the electronic repulsion 
is larger than the polaronic binding energy. For low vibrational modes 
the phase diagram is very rich. A noticeable point is the existence of 
a 4fc_F CDW phase for small values of the correlation strength. For 
realistic values of the electron repulsion and of the electron-phonons 
coupling constant a large phonons-mediated reduction of the Liittinger 
liquid parameter K p was found compared to the pure electronic model. 



I. INTRODUCTION 

Electron-phonons interactions have always attracted a lot of attention. Indeed, 
they are responsible for a large number of phase transitions, such as BCS supra- 
conductivity or Peierls transitions. Among the materials for which electron-phonons 
coupling are important, molecular materials constitute a special class. A trivial 
statement is that molecular crystals differ from simple crystals by the fact that 
their basic units are molecules. These basic units are therefore structured systems 
with a large number of internal degrees of freedom that may interact with the va- 
lence or conduction electrons. In particular, molecular crystals have two kinds of 
electron-phonons interactions. The first one is inter-molecular and the phonons 
essentially couple to the electronic structure through a modulation of the hopping 
parameter between two molecular sites. This electron-phonons coupling is responsi- 
ble for the Peierls transitions and has been extensively studied. The second type of 
interaction is intra-molecular. The totally-symmetric molecular vibrational modes 
couple to the electronic structure, essentially through a modification of the on-site 
parameters such as molecular orbital energies and on-site repulsions. In the 60's, 
Littlail suggested that intra-molecular vibrations could be responsible for supra- 
conductivity in organic conductors, more recently they were proposed as mediators 
for supra-conductivity properties in materials such as fulleridescl. 

Indeed in systems such as the organic conductors, there is aJwoad spectrum of 
intra-molecular vibrations ranging from 200cm~ 1 to 2000cm -1 Era, that is between 
~ t/5 and ~ t, where t is the hopping inter- molecular integral responsible for the 
conduction. Moreover it was shown that these vibrations couple quite well with 
the electronic degrees-jpf freedom both in the hight and low vibrational range of 
the Raman spectrumQu. This fact can easily be understood from simple chemical 
considerations. Indeed, the molecules acting as basic units in these systems share 
number of characteristics such are being large, strongly conjugated, and built from 
pentagonal cycles. These characteristics allow them to adjust their geometry to 
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their electronic charge at a low energetic cost, essentially by a modification of the 
angles in the pentagonal cycles. This mechanism leads to an important electron 
intra-molecular vibration coupling. 

Previous studies on correlated electronic systems coupled to intra-molecular vi- 
brations-have essentially explored the two asymptotic regimes, (i) the weak coupling 
regimeB'EI that can be treated by perturbative expansion from the electronic model 
and (ii) the strong coupling or polaronic regimeE3~t3. It is clear that the study 
of the intermediate coupling regime cannot be done by analytic treatments and 
requires up to date numerical techniques. The purpose of this work is to fill this 
gap and systematically explore the phase diagram of a one dimensional, correlated, 
electronic system coupled to intra-molecular vibrations, both as a function of the 
coupling constant and the correlation strength. The present paper systematically 
studies the Hubbard-Holstein model for a one dimensional quarter-filled chain for 
two values of the vibration frequency ,u> = 0.2i and lu = t, respectively correspond- 
ing to the top and the bottom part of the organic conductors Raman spectra. For 
each phase, properties such as spin and charge gaps, distance dependence of the 
spin, charge and singlet correlation functions, etc, are reported. 

The next section analyses the Hubbard-Holstein model and develops the compu- 
tational choices. The third section reports and discusses the results and the last 
section is devoted to the conclusion. 



II. MODEL ANALYSIS AND COMPUTATIONAL DETAILS 

A. The Hubbard-Holstein model 

The Hubbard-Holstein Hamiltonian associates a Hubbard Hamiltonian, which 
includes short range electron correlations, with dynamical phonons. The latter ace 
linearly coupled with the electronic degrees of freedom as in the Holstein modellij. 

H = H e + H p h + H e - V h 
with 

H e = e^Tli^ (4fl,cr t V + 4,o- C i+l,o-) + U ^ "i,T n i,J. 

i,cr i,a i 

H ph =cjJ2( b t b i + 1 / 2 ) 

i 

He-ph = .9 X! ni ( b i + 
i 

where c| CT , Ci il7 and the usual creation, annihilation and number operators 

of electrons of spin a on site i (iii — n^f -I- n^jj. b\ and &j are the intra-molecular 
phonons creation and annihilation operators. 

From the point of view of the isolated molecule the Hubbard-Holstein (HH) model 
tries to mimic the relaxation of the molecular geometry as a function of the ionicity. 
Indeed the on-site part of the HH model can be rewritten as 

Hi = en^ + U n l]nii + uj (b\bi + l/2 \ + gn z (b\ + b^j (1) 

(fit + m 9 -) ( bl + nA) +\] + (u- 2 £) n^i, + m (e - £) (2) 

The above formulation points out the three effects treated in the HH model. 
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1. The modification of the molecular orbital energy: e — ► e — g 2 /u>. This effect 
is very important (i) in multi-band systems since it strongly affects the relative 
filling in the different bands and (ii) in opened systems where electrons can 
jump in and out from an external bath. In our case it just changes the energy 
reference. 

2. The decrease of the effective on-site bi-electronic repulsion: U — ► U — 2g 2 /ui. 
In the strong coupling regime the effective electron-electron interaction can 
become attractive due to the electron-phonons interaction, and two electrons 
held together via molecular vibrations. 

3. The displacement of the harmonic oscillator describing the intra-molecular 



vibrations as a function of the molecular charge n*: 



1/2 



+ mg/u) (h + Hi g/uj) + 1/2 



This term mimics the relaxation of the 



molecular geometry as a function of the molecule ionicity. The A 



7- 



term acts as an effective molecular coordinate for which the equilibrium ge- 
ometry is linearly shifted from A = 0, when the molecular site does not carry 
any electron, up to A = 2g 2 /u>, when it carries 2 electrons. 

One sees from equation that the vibronic molecular states are coherent phonons 
states, eigenstates of the shifted harmonic oscillators. They can be referred as 
\rii, Szi, z^™ 1 ), where is the vibrational quantum number of the molecule i when 
it supports Hi electrons, 



\ni,Ssi,i/j*) (3) 



where 5 is the Dirac function. 



B. Computational details 

The calculations on the infinite chain are performed using the infinite system 
Density Matrix Renormalization Group methodta. The main problem risen by the 
HH model is the infinite number of vibronic states on each sites. In order to render 
the calculations feasible, the basis set have been truncated to the lowest vibronic 
states of each molecular site, that is \m, Szi, Vj) such that Vj = 0,1. This choice 
is physically reasonable since (i) we work at T = and therefore only the lowest 
vibronic states are expected to be involved, (ii) the molecules form well defined 
entities that are only perturbatively modified by the presence of their neighbors. 
The drawback of this choice is that we directly work in the vibronic basis set and the 
truncation of the basis set destroys any further possibility to separate the electronic 
degrees of freedom from the vibrational ones. We will see later, from a wave function 
analysis of the different phases, that the truncation does not affect the results of 
our calculations as long as we are not too close to a phase transition. 

In the phase diagram explorations we kept 100 states per renormalized block, 
while in the properties calculations 256 states were kept. The charge and spin gaps 
where computed using a double extrapolation (i) on the system size and (ii) on the 
number of states kept m. Typically the extrapolations over m were done from three 
DMRG calculations with respectively 100, 150 and 256 states kept. The maximum 
number of sites is 84 and the correlation functions presented below are done for 
this chain length. In some cases, exact diagonalizations on small clusters have also 
been performed in order to better analyze either the wave function or the energy 
spectrum. 
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III. RESULTS 



The present work explores in a systematic and unbiased way the whole range 
of electron-phonons coupling regime (from g/t — to g/t = 1.25) as well as the 
whole range of correlation strength (from U/t = to U/t = 16), for two values of 
the phonons frequency: u = 0.2t and to — t. These frequencies have been chosen 
in order to approximatively correspond to theJow and high part of the Raman 
spectrum of the TTFU, TMTSPa or M(dmit)j2 molecules. The adiabatic regime 
corresponds to the parameter domains u/t <C 1 and g/ui <C 1, and the strong 
coupling regime corresponds to g/w > 1 and u/t ^> 1. 

A. w = 0.2t 

Figure [l] reports the phase diagram for u = 0.2t as a function of g/u and U/t 
which are, with u/t, the three relevant parameters in the HH model (see eq. ||). 



Phase diagram o 1 2 
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g/ffl 

FIG. 1. Phase diagram of the Hubbard-Holstein model for ui = 0.2t. The 
solid line corresponds to U — 2g 2 /u> — 0, the long-dashed line corresponds to 
Ueff/t e ff = (U — 2g 2 /oj} jt exp(— g 2 /ui 2 ) = ±4, the dashed line corresponds to 
Ueff/teff = 600. The dotted-dashed liner-«orrespond to equal phonons polarization and 
electronic energies (as defined in referencalj). The inset blows up the small g/u and small 
U/t part of the graph. 

Five different phases have been characterized: two insulating phases where the 
electrons are strongly localized, a polaronic phase (diamonds) and a bi-polaronic 
phase (circles) and three metallic phases (crosses, stars and plus). These phases 
were characterized using the following set of tools 

(i) the variation of the energy with the number of sites, 

(ii) the spin and charge gaps, 

(iii) the spin, charge and singlet-singlet correlation functions, 
(iV) density matrices at the central sites. 
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1. The bi-polaronic phase 



In the strong coupling regime, U e ff = U — 2uj (g/ui) becomes negative and a bi- 
polaronic phase was found. As expected, the attractive character of U e ff strongly 
couples the electrons in pairs. For instance, for U/t — 4, g/t — 1 and 84 sites 
(U e ff = — 6t, g/u> — 5), the probability of having a lonely electron on the central 
site is smaller than 10~ 12 . The electron pairing induced by the intra-molecular 
vibrations is very strong, however, due to the Franck-Condon factors, this phase 
does not correspond to the singlet super-conducting phase but rather to a localized 
bi-polaronic phase. Let us suppose that an electron on a site i would like to hop 
on a neighboring j, omitting the spin degree of freedom the hopping term can be 
written as 



a\di |? 



Vg ) = t \rii - 1, v£ ; rij + 1, vl 3 ) 



7,<5 



n,-+l 



) \rk 



i 1 " 

Uj + l,v x 



(4) 
(5) 



that is the hopping integral between the molecular vibronic ground states is 
rescaled by the product of the Franck-Condom factors on the two sites i and j, 
( v a % \ v o { v o 3 \ v q 3+ )■ The relaxation energy or self-trapping energy (due to the 
vibrations) of the electron pair on a site can be evaluated as the difference between 
the vertical ionization potential (or electron affinity) and the adiabatic one, 



E 



relax 



(0 



For U/t — 4 and g/t = 1 this relaxation energy is as large as 5.00i. Figure ^ shows 
the overlap between the vibrational ground state corresponding to rij electrons on 
site i (v[] 1 ) and the different vibrational states corresponding to rij ± 1 electrons 
(i'2 i±1 )- One sees immediately that for large or even intermediate values of gju) the 
overlap between the low energy vibrational states is very small. The consequence is 
that the electron hopping between two neighboring sites is strongly hindered ; cither 
the transfer takes place towards small quantum-number vibrational states and the 
transfer integral is strongly renormalized by the small Franck-condom factors, or 
the transfer takes place towards large quantum-number vibrational states and it is 
hindered by the vibrational energetic cost. It is clear that the same phenomenon 
hinders even more the pair hopping since the displacement is twice as large and the 
Franck-Condon factor is squared ((v[] i \vq i±2 } = exp(— 2g 2 /u 2 )). 
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FIG. 2. Overlap between the ground-state of the undisplaced harmonic oscillator and 
the eigenstates of the displaced oscillator for different values of the displacement g/ui. 

If one looks at the energy per site as a function of the system size (or DMRG 
iteration number) , one sees that it is nearly constant. The amplitude of variation 
between the 4 sites system and the 84 sites system is smaller than 5 x 10~ 12 for 
U/t = 4 and g/u> = 5, the typical example we have chosen for this phase. This 
fact could have been surprising but it is easily understood in the context of a 
strongly hindered hopping. Taking into account the electron pairing and the strong 
localization one can reasonably approximate the system wave- function by tyg p 

TI TI TI 



of energy 

The quality of this approximation can be checked on the energy. The computed 
DMRG energy per site differs from E{^>bp) = —3.9 by at the most 7 x 1CT 12 in t 
units. Exact diagonalization on a 4-sites system, where the wave function can be 
explicitly analyzed, confirms this result. The projection of the exact ground state 
on ^ bp being larger than 1 — 10~ 7 . In fact the ground state of the 4-sites system is 4 
times quasi-degenerated and the infinite system ground state infinitely degenerated. 
Indeed, in the absence of inter-site Coulomb repulsion, any choice for the localization 
of the electron pairs is equivalent. Even with a 1/r Coulomb repulsion one sees that 
the infinite system should remain 4 times quasi-degenerated, the degeneracy lifting 
being of the order of magnitude of the rescaled hopping integral between the low 
energy vibrational states, that is of the order of t exp(— g 2 /ui 2 ) (1.4 x 10~ n for 
U/t = 4 and g/u> = 5). The charge gap is therefore exponentially small, the exact 
diagonalization of the 4-sites system shows that the first exited state with a "real" 
gap is a one-boson vibrational excited state at u). The spin gap is of the order 
of magnitude of U e ff (6.0 in our example) since it necessitates the breaking of an 
electron pair. 



2. The Liittinger liquid phase 

In the weak coupling regime — for small values of g/tu — up to the intermedi- 
ate coupling regime for intermediate values of the correlation-strength, one finds 
a phase which is essentially a Tomanaga-Liittinger liquidE-30, with parameters 
slightly rescaled by the presence of the vibrations compare to the purely electronic 
system. This result was expected, from continuity from the pure Hubbard model, 
and from previous works from Voit and SchulzQ (within a renormalization group 
(Rg) scheme in an incommensurate system). We will further refer to this phase 
as the Tomanaga-Liittinger (TL) phase. While restricted to a very small range 
of g/ui, for large values of the correlation strength, it is worth to notice that this 
phase extends up to values of g/uj larger than 1 for intermediate values of U (see 
figure [j]). For very small values of U, the competition between the TL phase and 
a bi-polaronic phase limits the extension of the former to a parameter range for 
which the effective repulsion remains positive, that is g/ui < v/Z7 /2u>. The study 
of the charge-charge and spin-spin correlation functions yield that the main effect 
of the vibrational degrees of freedom is to rescale the TL liquid parameters — in 
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particular the value of K p — compared to the purely electronic system. This result 
is in total agreement with the conclusions derived by Voit et aln. 

A simple approach would bring us to think that the value of K p is increased by 
the vibrational degrees of freedom since the effective electron-electron repulsion is 
strongly reduced by the electron-phonons interactions. Let us see what really comes 
out of the calculations. The correlation functions have been computed on a 84 sites 
system with 256 states kept in the DMRG procedure. They exhibit the expected 
behavior as a function of the inter-site distance with a power law decay, in agreement 
with the ungaped nature of both the spin and charge channel. Figure [| shows the 
charge structure factor for the pure Hubbard model as well as the Hubbard-Holstein 
one for the U/t = 1 and g/u> = 1.25 set of parameters. 

Charge structure factor 



m/t=0.2 - U/t=1 




FIG. 3. Charge structure factor for the pure Hubbard model (U/t = 1) and for the HH 
model (g/uJ = 1.25, U/t = 1). The small q part of the diagram is blown up in the inset. 



We can see on figure || that the numerical calculations contradict the simj 
prediction. The structure factor derivative at q = is directly proportional to K, 



2d 



d 



(6) 



Kp= fq\ Jo <(ne(0) " (ne(0)>) (ne(r) ~ {ne{r))))eiqT 

and is clearly smaller for the U/t = 1 and g/u) — 1.25 set of parameters than for the 
U/t = 1 Hubbard model. The K p values computed from the structure factors are 
0.79 for g/uj — 1.25 and 0.89 for the Hubbard Hamiltonian — in total agreement 
with the-Kalues found in the literature from numerical resolution of the Bethe ansatz 
solutionis. This unexpected reduction of K p was already noticed by Voit et aln 
from Rg considerations. One should notice that it cannot be attributed to the 
rescaling of the hopping integral since 

- - ^ CD 

ttff texp(-g z /uj z ) 
and the corresponding value of K p is much larger than 0.850. 



3. The 4k F CD W phase 

In a regime of intermediate coupling and intermediate to large correlation 
strength, we found a phase for which the charge correlation functions are domi- 



7 



nated by 4kp charge density wave (CDW) fluctuations. The multiplication of the 
correlation functions by r 2 (r being the inter-sites distance) allows the elimination, 
when the Fourier transform is performed, of the 1/r 2 term that contribute at all 
frequencies. The frequency analysis of r 2 xthe correlation function is therefore much 
clearer, and the relative importance of the 2kp and Akp terms enhanced. Figure ^| 
reports this Fourier analysis for U/t = 4 and g/uo — 1.25. One can see that while 
the spin-spin correlation function is change very little compared to the Hubbard 
model, the 2kp contribution to the charge-charge correlation function has disap- 
peared and a strong 4kp contribution has set place. Figure [| reports the same 
correlation functions as a function of the inter-sites distance. Both charge and spin 
correlation functions decrease as a power law as a function of the inter-sites dis- 
tance speaking in favor of a gap-less system in both the spin and charge channels. 
The charge and spin gaps have been computed independently from double extrap- 
olations (i) on the chain length and on the number of states kept in the DMRG 
calculations, i.e. 100, 150 and 256. While the charge channel clearly extrapolates 
toward a null gap, in the spin channel the question is not as clear. Indeed, the 
gap is found to be A CT (100) = 8.1 x 10~ 3 for a DMRG calculation where 100 states 
are kept, A CT (150) = 8.3 x 10~ 3 for 150 states kept and A CT (256) = 8.6 x 10~ 3 for 
256. These gap values are very small however they do not extrapolate toward a 
null value when the quality of the calculation increases. This behavior pleads for a 
non null but very small gap in the spin channel. One should note that this result is 
not incoherent with the spin-spin correlation function behavior since such a small 
gap means a very large coherence length, of the order of magnitude of A^ 1 , i.e. 
the exponential behavior of the correlation function should take place at inter-sites 
distances larger than the chain length. 
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FIG. 4. Fourier transform of r 2 
times the charge-charge (b) and 
spin-spin (a) correlation functions 
for the set of parameters U/t — 4 
and g/u = 1.25. The angle is given 
in degrees. 
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FIG. 5. Absolute value of 
charge and spin correlations 
as a function of the inverse 
of the inter-sites distance 
(U/t = 4 and g/u = 1.25). 



In the Luttinger liquid theory the Akp phase is supposed to occur for values of the 
K p parameter smaller than 1/3, that is (i) in very strongly correlated systems and 
(ii) for values of K p unreachable in the Hubbard model — for which 0.5 < K p < 1. 
The presence of intra-molecular vibrations not only allows the existence of a Akp 
CDW phase, but this phase can be reached for values of the bi-electronic repulsion 
as small as U/t = 2.5 (and most probably even lower). We computed the K p values 
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from the structure factors and found that similarly to what happened in the TL 
phase, the K p value is strongly diminished in comparison to the purely electronic 
case. For U/t = 4 and g/t — 0.25 the value of K p is 0.56 while it is 0.71 for the 
U/t = 4 Hubbard model. One should however note that the 0.56 value is in good 
agreement with the K p expected for the effective parameters, U e ff/t e ff = 16.1. 

4- The small polarons phase 

For positive U e ff but relatively large values of g/u), the system is in a small 
polarons phase (diamonds on fig. |l|). Indeed the electrons are no longer paired 
but are still strongly localized. The probability of having two electrons on the 
central site is smaller than 10 -9 for U = 8, g/t = 0.75 and 82 sites (g/u> = 3.75 
and U e ff = 2.375). These lonely electrons remain strongly coupled to the intra- 
molecular vibrations. The on-site vibrational relaxation energy is as large as 2. St 
and the rescaled hopping integral between the ground vibrational states as low as 
7.8 x 10~ 7 . The ground-state wave-function can be approximated by the totally 
localized wave-function \Pp 

T i T I T 



and the residual hopping term treated perturbatively. The overlap of the exact 
ground-state wave-function of an 4-sites system with "fp is as large as 0.99 for 
U/t — 8 and g/t = 0.75 and even as large as 0.98 for U = 4 and g/t — 0.5 which 
is very closed to the phase transition. Coherently, the DMRG energies per site 
remains nearly constant as the system size increases, with a maximal variation of 
3 x 10~ 7 from E(^ P )/N sUes = -l/2g 2 /uo + l/2u), for U = 8 and g/t = 0.75. 
As for the bi-polaronic phase the absence of Coulombic inter-site repulsion in our 
model induces a strong quasi-degeneracy of the ground-state: 12 times for the 4- 
sites system but infinitely for the infinite system. Even in the presence of a 1/r 
repulsion, the system would remain at least twice quasi-degenerated due to the 
even sites versus odd sites equivalence with an additional factor of 4 due to the spin 
quasi-degeneracy. Capone et alt2 have studied the conditions of existence of this 
phase by exact diagonalization on small clusters (4 and 6 sites). They found that 
the localized polarons phase is stable as long as the polarization energy defined as 
£poi = —N S it es p(l — p)g 2 /&> is larger than the electronic energy e e ; ec . p is the filling 
of the system. The curve e po i / e e iec = 1 is reported on fig. [j] with a dashed-dotted 
line and fits relatively well with the phase boundary in the infinite system. In the 
same spirit curves of constant U e ff/t e ff have been plotted. This small polarons 
phase corresponds to very large values of the f7 e // /teff ratio. 

5. The Luther-Emery phase 

Finally for negative values of t/ e //, but large values of t e ff a delocalized phase 
takes place for which all spin, charge and on-site singlet correlation functions seem 
to decrease with the inter-site distance as a power law. Figure || reports the absolute 
values of the correlation functions for U/t — 0.05 and g/t = 0.1, that is g/uo = 0.5, 
U e ff = —0.05 and t e ff — 0.78. In all computed cases (U/t — 0.2, g/u> — 1 and 
1.3 ; U/t = 0.05, g/oj = 0.5) the charge density fluctuations dominate, however as 
Ueff/teff decreases the on-site singlet fluctuations increases compared to the charge 
ones. A direct calculation of the charge and spin gaps yield a clearly ungaped 
charge channel and a slightly gaped spin channel. The computed spin gaps are 
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A CT = 0.05 for U/t = 0.05 and g/t = 0.1, and A CT = 0.06 for U/t = 0.2 and 
g/t = 0.2. These very small values are compatible with the the power law decrease 
of the spin correlation functions since with such small gaps, the correlation length 
is very large — of the order of magnitude of A^ 1 — and therefore the exponential 
decay of the correlation function cannot take place at distances smaller than the 
order of magnitude of the chain length. One can therefore identify this phase with 
a Luther- Emery modelEI The values of K p extracted from the structure factors are 
again smaller than the values of the purely electronic model, with K p = 0.99 for 
U/t = 0.05 and g/t = 0.1, and K p = 0.89 for U/t = 0.2 and g/t = 0.2. Of course, 
they are as well smaller than the K p values corresponding to the negative U e ff . 

Correlation functions 

<a/t=0.2 U/t=0.05 g/t=0.1 




FIG. 6. Absolute value of the charge (stars), spin (crosses) and singlet (plus) correlation 
functions for U/t = 0/05, g/t = 0.14 and u/t = 0.2 



B. u = t 

Figure [?] reports the phase diagram for u — t as a function of g/uj and U/t. 
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FIG. 7. Phase diagram of the Hubbard-Holstein model for u = t. The solid line corre- 
sponds to U — 2g 2 /u — 0. 



For this intermediate value of the phonons frequency, one finds only three phases, 
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a slightly perturbed Liittinger liquid phase for U e ff > 0, a Luther-Emery phase for 
U e ff < and small g/uj, a localized bi-polarons phase for U e ff < and larger g/u). 

For U e ff > 0, the phase is delocalized, dominated by 2kp SDW fluctuations. 
Both charge and spin correlation functions decrease as power laws as a function 
of the inter-sites distance and both spin and charge channels are ungaped to nu- 
merical accuracy. The structure factors for the U/t = 4 Hubbard model and for 
the Hubbard-Holstein model with U/t — 4 g/u = 0.5 and U/t = 4 g/u) = 1 arc 
reported in figure ||. As can be seen the HH model is in this case indistinguishable 
from the pure Hubbard model and, as a consequence, the K p values as not sensitive 
(to numerical accuracy) to the electron-phonons coupling. In fact everything goes 
as if we were in the adiabatic regime while the model parameters can be as far from 
it as u>/t — 1 and g/u> — 1.25. This is typically a case where the Born-Oppenheimer 
approximation is fully valid and the electronic and vibrational degrees of freedom 
are nearly independent. 
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FIG. 8. Structure factors for U/t = 4. The solid line corresponds to the pure Hub- 
bard model, the crosses and the circles to the Hubbard-Holstein model with respectively 
g/t = .0.5 and 1. 

For U e ff < one has a Luther-Emery phase for small values of g/cu. This phase 
is very similar to the one found for u/t = 0.2, with an effective attractive on-site 
interaction and large effective hopping integrals. The CDW and the on-site singlet 
fluctuations dominate, and become of the same order of magnitude for g/io = 0.5 
and U/t = 0.2. For larger values of the electron-phonons coupling, the CDW are 
larger than the singlet fluctuations (g/u) — 1 and U/t = 1). As for the u> = 0.2 case, 
all three correlation functions — spin-spin, charge-charge and singlet-singlet - 
present a power law behavior as a function of the inter-sites distance for the 84 
sites system. The direct calculation of the charge and spin gaps yield and ungaped 
charge channel and a slightly gaped spin channel with A CT = 0.04t for U/t = 1 
and g/t — 1, as well as for U/t — 0.2 and g/ui = 0.5. Again, the smallness of the 
spin gap explains the power law behavior of the spin-spin correlation functions at 
the computed distances. This time however the values of K p are larger than in 
the Hubbard model, in agreement with the attractive effective on-site interaction. 
K p ~ 1.1 for both U/t = \g/t= \ and U/t = 0.2 g/t = 0.5. 

When g/uj increases — but still for negative U e ff — one goes from the Luther- 
Emery phase toward a localized bi-polaronic phase through a cross-over. Indeed, 
for U/t — 1 and g/u = 1.5 the system is clearly localized in a bi-polaronic phase 
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with the usual characteristics : nearly constant energy per site, large projection of 
the wave function on f bp (0.94 for the 4 sites system), etc... For U/t = 1 and 
g/oj = 1.25 however, the system is in an intermediate regime. The derealization is 
still reasonably large with an effective hopping of 0.21i and a projection of the 4 sites 
system wave function on ^bp of only 0.72. In this cross over regime, the system is 
gaped both in the spin and charge channels, with a spin gap of A CT = 0.43 and a 
charge gap of A p = 0.84. Figure ^ reports the charge, spin, and singlet fluctuations. 
As expected all of them decrease exponentially with the inter-site distance. When 
the electron-phonons coupling increases, the system localization increases and the 
wave function tends toward the totally localized wave function^sp. Coherently, 
the spin gap increases and is expected to follow the same variations as |£7 e //| = 
2g 2 /us — U, since one needs \U e ff\ in order to break the singlet pairing on the 
sites. In totally localized systems (g = +oo), the charge channel is ungaped since 
in the absence of any derealization integral, the different possible choices for the 
localization of the pairs have the same energy. The charge gap can therefore be 
expected to be scaling as the rescaled hopping integral t e ff = texp(—g 2 /io 2 ). That 
is exponentially decreasing when the electron-phonons coupling increases. 

Correlation functions 
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FIG. 9. Absolute values of the charge (stars), spin (cross) and singlet (plus) correlation 
functions in a semi-log scale. U/t = 1, g/t — 1.25 and id = t. 



IV. CONCLUSION 

We have investigated the phase diagram of the Hubbard-Holstein Hamiltonian in 
a quarter-filled chain for two values of the phonons frequency. These frequencies 
have be chosen in order to correspond to the low (lo — 0.2i) and hight (lo = t) 
frequencies of the intra-molecular totally-symmetric vibrations of TTF, TMTSF 
and related molecules from which the organic conductors are built. 

For the lo = t frequency, the phase diagram is very simple. For a positive effective 
on-site interactions {U e ff = U — 2g 2 /uj > 0), the system is in a Liittinger-Liquid 
phase with slightly rescaled parameters compared to the purely electronic system. 
For negative effective on-site interactions (U e ff < 0) and small electron-phonons 
coupling, the system is in a Luther-Emery phase, however the spin gap remains 
very small and the coherence length large. The K p parameter is again rescaled 
compared to the purely electronic system, in agreement with the effective model 
U e ff — U — 2g 2 /uj, t e ff — texp(~g 2 /uj 2 ). When the electron-phonons coupling 
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increases, the system goes toward a localized bi-polaronic phase through a soft 
cross-over. 

At low frequency (uj — 0.2t) however, a rich phase diagram have been found with 5 
different phases, some delocalized and some strongly localized. For positive effective 
interactions one has three different phases. For small electron-phonons coupling and 
intermediate coupling for small correlation strength, the system is in a Luttinger 
Liquid phase. For intermediate electron-phonons coupling, the system is in a phase 
for which the 2kp charge fluctuations disappear and the charge correlation functions 
have only 4/cf CDW fluctuations. For large electron-phonons coupling, the system 
is strongly localized in a polaronic phase. When the effective interaction become 
negative, one finds a Luther-Emery phase for small electron-phonons coupling, and 
a strongly localized bi-polaronic phase for intermediate to large couplings. The 
Luther-Emery phase has a very small spin gap (of the order of 0.05i and the spin- 
spin correlation functions exhibit a power law decrease up to at least 60 to 80 sites) . 
Unlike for w/t = 1, for u>/t = 0.2 the delocalized phases do not behave as would the 
U e ff = U — 2g 2 /uj, t e ff — t exp(— g 2 /ui 2 ) effective Hubbard model. Indeed the K p 
parameter is strongly decreased due to the electron-phonons interactions while the 
decrease of the effective repulsion (compared to the pure Hubbard model) would 
have led us to expect an opposite behavior. 

This surprising decrease of the K p parameters should be put in context with val- 
ues of the K p parameters found in photo-emission experiments on TMTSF and 
TMTTF quasi-one-dimensional organic compoundsEirEl Indeed Zwick et alt3 
found a density of states exponent a — 1/4 (if p + 1/K p — 2) slightly larger than 
1 (K p < 0.25) in the metallic phase. According to the Luttinger Liquid theory the 
system should be insulating fox. such large values of a. Tentative explanations have 
been made by Schwartz et at3, involving an effective doping of a Mott insulator 
due to the inter-chain hopping. On another hand, the electron intra-molecular vi- 
brations interactions are never considered in these systems. For the low frequency 
modes of these systems, the typical energy scale of the electron-phonons interac- 
tion : 2g 2 /u! is of the order of magnitude of one-half to one-fourth the intra-chain 
hopping amplitudeEm, that is in the intermediate regime of our phase diagram. For 
such values of the electron-phonons coupling constants and values of the electronic 
repulsion between one to four t, the system is in a Luttinger liquid phase for which 
the intra-molecular vibrations strongly decrease the K p parameter compared to its 
electronic value. It is clear that the Hubbard model lacks at least nearest-neighbor 
bi-electronic repulsion in order to well represent these systems. It can however be 
expected that the trend toward a strong reduction of the K p parameters, observed 
in the Hubbard-Holstcin model, would be similar with longer range interactions, 
and that the coupling of the electronic degrees of freedom with the intra-molecular 
vibrations would lead to values of the a parameter unreachable in a purely elec- 
tronic model. Considering these trends, it is reasonable to think the the electron 
intra-molecular vibrations coupling will have to be considered in order to obtain a 
realistic description of the observed photo-emission experimental data. 
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